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Abstract During this work, an interacting chameleon like 
scalar field scenario, by considering SNela, CMB, BAO and 
OHD data sets is investigated. In fact, the investigation is re¬ 
alized by introducing an ansatz for the effective dark energy 
equation of state, which mimics the behaviour of chameleon 
like models. Based on this assumption, some cosmological 
parameters including Hubble, deceleration and coincidence 
parameters in such mechanism are analysed. It is realized 
that, to estimate the free parameters of a theoretical model, 
by regarding the systematic errors it better the whole of the 
above observational data sets to be considered. In fact, if 
one considers SNela, CMB and BAO but disregards OHD 
it maybe leads to different results. Also to get a better over¬ 
lap between the counters with the constraint < l.thexl 
function could be re-weighted. The relative probability func¬ 
tions are plotted for marginalized likelihood ^ (i2nio, Oli, j3) 
according to two dimensional confidence levels 68.3%, 90% 
and 95.4%. Meanwhile, the value of free parameters which 
maximize the marginalized likelihoods using above confi¬ 
dence levels are obtained. In addition, based on these calcu¬ 
lations the minimum value of based on free parameters 
of an ansatz for the effective dark energy equation of state 
are achieved. 

PACS 98.80.-k-98.80.Es-95.36.H-x 

1 Introduction 

Observational data sets including Cosmic Microwave Back- 
ground(CMB) [1,2], Supernovae type la (SNela) [3,4], Bary- 
onic Acoustic Oscillations (BAO) [5,6], Observational Hub¬ 
ble Data (OHD) [7,8], Sloan Digital Sky Survey (SDSS) [9, 
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10], and Wilkinson Microwave Anisotropy Probe (WMAP) 
[11,12], are considered as criterion for accuracy of theo¬ 
retical models. Amongst these constraints, the CMB and 
SNeIa( because of abundance of their data sources) attract 
more attention. It is notable, the SNela constraint in high 
redshift values do not give a good clue to investigate the 
evolution of the Universe. It is obvious the results of in¬ 
dividual observations give different values for free param¬ 
eters of a theoretical model; hence, it is better that, to es¬ 
timate the best quantities for free parameters of the model 
one considers the whole of observational data sets including 
CMB, SNela, BAO and OHD. Therefore, this motivated us 
to study the behaviour of free parameters and their overlaps. 
Thence, a collective of observations including SNela, CMB, 
BAO and OHD are considered. Meanwhile the mentioned 
observational data sets have predicted an ambiguous form 
of matter which leads to an accelerated phase of present 
epoch and it is well known as dark energy. Based on this 
ambiguous form of matter, scientists have proposed differ¬ 
ent proposals up to now. Amongst all of those proposals, 
the cosmological constant. A, model attracts more attention 
[13,14]. But this mechanism suffers two well known draw¬ 
backs. The first of them is related to estimate the contribu¬ 
tion of quantum fluctuation of zero point energy and second 
is related to the ratio of A and dark matter energy densi¬ 
ties. These problems and also the excellent work by Brans 
and Dicke [15] motivated scientists to introduce a mecha¬ 
nism which A had time dependency, namely quintessence 
[16, 17, 18]. Beside the quintessence mechanism, some pro¬ 
posals which risen from quantum gravity or string theory 
are introduced to estimate cosmological parameters. For in¬ 
stance one can refer to tachyon [19,20], phantom [21,22, 
23], quintum [24,25], k-essense [26,27] and so on. Also 
some models which were risen from quantum field fluctu¬ 
ations or space time fluctuations attract more attention to 
investigate dark energy concept. For such models, one can 
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mention Zero Point Quantum Fluctuations(ZPQF) [28,29, 
30], Holographic Dark Energy(HDE) [31,32,33,34,35,36], 
Agegraphic Dark Energy (ADE) and new-ADE [37,38,39]. 
If scalar field, in the quintessence model couples to mat¬ 
ter (non relativistic) it induces to appear a fifth force. When 
the coupling is of order unity, the results of strongly cou¬ 
pling scalar field is not in good agreement with local grav¬ 
ity tests (for instance in solar system). Thus a mechanism 
should be exist to suppress the effect of fifth force; such 
mechanism is capable to reconcile strong coupling models 
with local experiments was proposed by Khoury and Welt- 
man [40,41] and also, separately, by Mota and Barrow [42] 
namely chameleon. In this mechanism, one can not choose 
an arbitrary Lagrangian for matter, L,„. To avoid deviation 
of geodesic trajectory, the author of [43] has shown the best 
choices are L,„ = P and L,„ = —p, where P is pressure and 
p is energy density of matter, for more discussion we re¬ 
fer reader to [44,45,46]. Therefore main motivation of this 
work is investigation the behavior of an interacting scalar 
field mechanism; based on these calculations and SNela, 
CMB, BAO and OHD data sets the minimum value of 
for the effective dark energy equation of state are achieved. 
The scheme of the paper is as follows: 

The above discussions which are brief review about obser¬ 
vational and theoretical motivations are considered as intro¬ 
duction. In Sec. 2, the general theoretical discussions risen 
from a chameleon like mechanism related to the cosmologi¬ 
cal parameters such as Hubble, deceleration and coincidence 
parameters will be discussed. In Sec. 3, a brief review about 
cosmological data sets are brought. In Sec. 4, the observa¬ 
tional data sets including SNela, CMB, BAO and OHD are 
considered, to estimate the minimum value of x'^ related to 
free parameters of an ansatz for the effective dark energy 
equation of state. And at last. Sec. 5, is dedicated to con¬ 
cluding remarks. 


2 Conservation and field’s equations in an effective 
dark energy scenario 


and matter sector. It should be noted L is the Lagrangian den¬ 
sity of matter which consists of both dark matter and dark 
energy sectors as perfect fluid [47,48,49,52,53]. It should 
be noticed that, background is a spatially flat Eriedmann- 
Limatire-Robertson-Walker (ELRW) Universe, with signa¬ 
ture (-1-2). The variation of the action 1, with respect to 
(w.r.t) guv results the gravitational field equation as 

Gmv=/(<P)7mv + 7'^v\ (2) 


where the stress-energy density of scalar field expresses 
= (^'^^l(P^v(p-^8^lv{^(pf^-g^vV{(p), (3) 


and 


- ^ -2 5{vf^L) 


(4) 


is the definition of stress-energy tensor of matter. By consid¬ 
ering 00 and ii components of the energy density and 
pressure could be achieved. After some algebra the conser¬ 
vation equation reads 






= 0 . 


(5) 


In addition, the variation of the action 1, w.r.t scalar field 
gives the evolution equation as 


(p + 3H^ 


-V{(p) + 


dm, 

d(p 


(6) 


Now, by substituting equation 6 into relation 5, two conser¬ 
vation equations for scalar field and matter are attained as 


v^[/((p)r^0] = -/((p)L, 


(7) 

(8) 


where over dot denotes derivation w.r.t ordinary cosmic time, 
f. As it was mentioned in the introduction, the Lagrangian of 
matter is considered as L = -f , [52,54] ; where sub¬ 
script m denotes matter (cold dark matter and baryons) and 
de refers to dark energy. Then the conservation equations 
could be rewritten as 


In chameleon like scalar field scenario, the mass of scalar 
field is a function of local matter density, so that, it is suf¬ 
ficiently large on dense environment. Due to this fact, the 
equivalence principal (EP) is satisfied in the laboratory [40, 
41]. In addition, the Brans Dicke (O parameter for two ob¬ 
servational values of 7 post Newtonian parameter take the 
values of order 10“* [47], which satisfies the solar system 
constraint. The chameleon like scenario is defined as 

5 = J d^x^m8(R-^^^(p^^(p-2Vi(p)+2f{(p)Ly (1) 

[43,44,48,49,50,51,52]. In this equation g is the determi¬ 
nant of the metric, V (^) is a run away potential and the latest 
term indicates a non-minimal coupling between scalar field 





(9) 


nmS) 

= -f{<P)L(de)^ 

( 10 ) 


V 

= f((P){L(m) +i-(de))- 

( 11 ) 


By combining the relation 2 and above equations, it is easy 
to receive 


\<P)' 


/(<P)7'(S 




( 12 ) 


In the next step, by virtue of the definition of and equa¬ 
tion 2, the Einstein tensor is modified as 


=/(<?) 


^ (m) 




gv 

(de) 


_J_jgV 

fi(p) 


(13) 
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Hereafter, we postulate that both scalar field and dark en¬ 
ergy, behave the same as perfect fluid, thence for such per¬ 
fect mixture the effective stress-energy tensor is obtained as 
follows 

where subscript DE denotes effective dark energy. There¬ 
fore using equations 7-14, the modified Einstein equation 
and conservation relations are attained as 




rpliV rp^V 




m)T, 


(m) 

MO 

(DE) 


= -f{(p)L(,n), 


(15) 

(16) 
(17) 


It should be noticed that, in the right hand side of above 
equations, only is appeared. In fact it could be con¬ 
cluded that the energy, for different components of the Uni¬ 
verse is not conserved separately. In refs. [55, 56, 57], it has 
been shown that, for perfect fluids, that do not couple di¬ 
rectly to the other components of the Universe, there are 
different Lagrangian densities are equivalent. Namely, one 
can find that the two Lagrangian densities = P and 
L{m) = ~P give the same stress-energy tensor and the equa¬ 
tion of motions for all components of the system are simi¬ 
lar as well. But in an interacting case, which matter has an 
interaction with scalar field, the Lagrangian degeneracy is 
broken. Based on ref. [43], the best choice for such models 
is = P. Using this definition for Lagrangian of the mat¬ 
ter one can obtain 


= \f{^) 

and also 


Pm + PdE 


fi(p)Pm +3Hf{(p)p,n = 0, 


f{(p)pDE +3Hf{(p) 


Pde = 0, 


(18) 

(19) 

( 20 ) 


1 + (Ode 

where H = d{t)/a{t) is the Hubble parameter, a{t) is scale 
factor and (Ode, is the EoS parameter of the effective dark 
energy and satisfies EoS equation as 

Pde = (Ode x Pde- (21) 

To establish an accurate link between theoretical results and 
observations, one can use the red shift parameter, z, instead 
of the scale factor; these two cosmological parameters have 
a relation as 
a (to) 


i{t) 


= 1 +z z = —(1 +z)H. 


( 22 ) 


Thus substituting equation 22 into 19 and 20, one finds out 
./(^)Pm=/oxp,„ox(l-b^)^ (23) 


f{(?)PDE = /o X Pdeo X exp [ 


(24) 


where Pdeo and p„,o refer to energy densities of dark energy 
and matter at present time, respectively. 


2.1 Hubble parameter 

Dimensionless Hubble parameter and density parameters could 
be defined as 


Hq 

(25) 

PmO 

37/2’ 

(26) 

/s Pdeo 

- 37/2 ■ 

(27) 


The dimensionless density parameters could be rewritten as 

= (28) 

5 Hq 

iiDE0 = ^^- (29) 


foXpDEO 

Therefore using relations 18 and 25, the dimensionless Hub¬ 
ble parameter is obtained as follows 

3i±^430) 


e '^{ z ) = A«o(l +z)^ + f2D£oexp 




1 +z 


2.2 Coincidence parameter 


The ratio of dark matter and dark energy is defined as coin 
cidence parameter and could be obtained as 

^ _ Pm 

Pde 

1 + (Ode{z) 


(31) 


I n 3 r o r( + (ODE[Z) -- 

(l-fz)^exp-3 / ———— dz]. 
Jo 1 +z 


(z)- 


(32) 


= fo 

Also one can obtain 
^ _ -3code{z) 
dz 1 +z 

Due to the role of this parameter, r, in the investigation of the 
cosmic evolution, it attracts more attention in observational 
investigations. In fact one can observe that, this importance 
is arisen from the relation between the EoS parameter and 
the evolution of r. 


2.3 Deceleration Parameter 

To investigate the acceleration of the Universe, one can use 
deceleration parameter which is defined as 

The above equation can be rewritten as 

3 f {1 +COde)E^-Q mo{^+z)^0)DE\ 

= + - )■ O'** 

In present epoch of the Universe evolution, deceleration pa¬ 
rameter is determined as 
1 3 r. 


— 2 2 ^ EimO 


(Ode (0). 


(35) 
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To solve the above equation we introduce an ansatz for EoS 
parameter as [54,58] 

®D£(z) =-1+®o+®i(l+2 )^, (36) 

where coq, (0i and p are free parameters of the model, where 
the minimum value of of them will be obtained in fit¬ 
ting part. Also it is notable if we choose P —0, the model 
reduces to EoS constant models (for instance ACDM) [54]. 
By substituting 36 in 30, the dimensionless Hubble parame¬ 
ter is attained as follows 


E^[{z\Vi})=a,nQ{\+zf + aDE0{\+zf^ X 
exp 


(37) 


^3^ 


where 




(38) 


and {P;} is a set of free parameters which should be deter¬ 
mined using data fitting process. Using equation (36), one 
can rewrite the equations 23, 24 and 31, respectively as 


f{(p)pm 

fi(p)pDE 


foPmO X (1 

foPOEO X (1 X 


exp 




(39) 

(40) 


and 

riz) = ro(l +z)^^^^^“°^exp 




(41) 


3 A brief review as to cosmological observational data 
sets 

In this section, we should emphasis that the analysis is re¬ 
stricted to the background level, and do not include pertur¬ 
bations. In the following, we want to compare our theoreti¬ 
cal results with observations. To this end, we consider four 
important data sets including SNela, CMB, BAO and OHD. 
In some papers, it was claimed OHD, which obtained versus 
red shift, is comparable with SNela data set, for instance we 
refer reader to reference [7] and references which are there. 
This subject motivated us to investigate the effects of this 
new data set beside other observations to improve the the¬ 
oretical results. As it will be discussed, the results of OHD 
although is not independent of SNela and BAO data sets [7] 
but has not any dependency to CMB. Also there are two 
ways to study CMB and BAO data point among the full pa¬ 
rameter distribution and Gaussian which in follow the latter 
will be used. 


3.1 Supernovae type la 

It is explicit that, supernovae attract more attention in em¬ 
pirical cosmology. Whereas they are very luminous, people 


interested to consider them, also for instance at closer dis¬ 
tances (i.e. lower redshift) they could be used to calculate 
Hubble parameter, and for farther distances (i.e. higher red- 
shift) they attain an important role to estimate deceleration 
parameter q. It is obvious there are uncertainties of different 
nature; statistical or random errors and systematic errors. In 
this work it is remarkable the systematic errors for SNela 
and OHD are neglected. In reality there is always a limit 
on statistical accuracy, besides the trivial one that time for 
repetitions is limited. The assumption of independence is 
violated in a very specific way by so-called systematic er¬ 
rors which appear in any realistic experiment. Eor instance 
experiments in nuclear and particle physics usually extract 
the information from a statistical data sample. The precision 
of the results then is mainly determined by the number N 
of collected reactions. Besides the corresponding well de¬ 
fined statistical errors, nearly every measurement is subject 
to further uncertainties, the systematic errors, typically as¬ 
sociated with auxiliary parameters related to the measuring 
apparatus, or with model assumptions. The result is typically 
presented in the form 

X = 2.34±0.06 = 2.34±0.05 {stat) ± 0.03{syst). 


The only reason for the separate quotation of the two uncer¬ 
tainties is that the size of the systematic uncertainties is less 
well known than that of the purely statistical error [60]. By 
virtue of the likelihood functions, one able to estimate the 
minimum value of for the set of parameters {pi}, as 


^({Pi,Mo}) °=exp 


' 1 7 

-^XSNeiiPhPo}) 


(42) 


where 


/■(„ ,,r\\ [P(ibs{Zn) — PlhiZn', {Pn Po})] 

XsA'edPii Mol) 2 ^ 2 

n=\ 


(43) 


In 43, Pobs{zn) is the observational distance modulus for 
nth supernova, a„ is the variance of the measurement and 
Pth{zn) is the theoretical distance modulus for nth supernova 
which defined as 


Mr/.(2«;{Pi,Mo}) =51ogio[DL(Zn;{pi})]+)iO, 
^0 =42.38-51ogio [h], 

dz 


£>L(z«;{Pi}) = (1 


d / 


E(z;{pi})’ 


where Dl is the luminous distance and h = l OOkms^ ^Mpc^ *. 
To achieve best fit of free parameters, one can marginalize 
likelihood function w.r.t /io [59,60]. Thence x'sNeiiP'}) 
duces to 

X|ve({Pi}) ~ (44) 

where A, B and C are defined as follows 

^ ^ [Pohs jZn ) ~ Pth (Zn i (Pi! Mo ~ 0} )] 
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n l^obs(Zn) ^J^th{Zn7{Pi^ — 0 }) 

B=L - ^2 -> 

n=l 
557 j 


^=L- 2 - 


(46) 

(47) 


3.2 Cosmic Microwave Background 


According to oscillations appear in matter and radiation fields 
Doppler peaks in radiation (photon) spectrum are produced. 
Also it should be noted that existence of dark energy, affects 
the place of the Doppler peaks in spectrum diagrams. To 
determine the shift of these peaks, theoretically, CMB shift 
parameters is defined as refs. [1,61] 


{Zrec '7 {Pi}) 


^mO dz 

fo Jo £(z;{pi})' 


(48) 


In CMB investigations [62], the Xcmb function versus CMB 
shift parameter is 


Xcmb{{V'i\) 


[f^obs {Pi})] 


(49) 


where Rghs = 1.725, <Jr = 0.018 and Zrec ~ 1091.3 are ob¬ 
servational quantities of CMB shift parameter, uncertainty 
of R in C7i confidence level and recombination redshift, re¬ 
spectively refs. [61,1]. 


3.3 Baryonic Acoustic Oscillations 


As in [63] mentioned, because BAO can be considered as 
a standard length scale at a wide range of redshift it is an 
useful candidate for cosmological models testing. The im¬ 
portance of BAO mechanism is related to its ability in esti¬ 
mation the contents and curvature of the Universe. One can 
establish a relation between theoretical BAO parameter, A,/„ 
and dimensionless Hubble parameter, Eq.(30), as 


A,h{zb-,{Pi}) 


(50) 



1/3 


J _ rZb dz 
Zb -lo £(f;{pi}) 



where Zb = 0.35 [5,6]. Also Xbao BAO mechanism in¬ 
vestigation is as follows 

= (51) 

and also Aobs = 0.469(n,/0.98)-‘’'35 and = 0.968, [6,59]. 
It is obvious that, the BAO are detected in the clustering 
of the combined 2dFGRS and SDSS main galaxy samples, 
and measure the distance-redshift relation at z = 0.2. But we 
consider BAO in the clustering of the SDSS luminous red 
galaxies in which measure the distance-redshift relation at 
z = 0.35 [64]. 


3.4 Observational Hubble Data 


We suggest that, if people want to investigate the accuracy of 
any theoretical model, it is better, maybe, to consider SNela, 
CMB, BAO and OHD together. In [7], it was claimed that 
three different models of dark energy i.e. ACDM, ^CDM 
and XCDM have been investigated just by considering H{z) 
measurement. But they have used Hq = 68 ± 2.8 and Hq = 
73.8 ± 2.4 which risen from SNela data [8]. Therefore it is 
realized that for the comparison between theoretical results 
and observations only OHD could not be considered. The 
XoHD function parameter based on OHD data set is defined 
as 


XoHD({Pi>77o}) 


^ [Hobs{Zn)-HoEth{z„-,{Pi})? 


(52) 


after marginalize w.r.t Hq, to calculate likelihood, Xohd could 
be considered as 

ZoHD({Pi}) = A//— (53) 

where 


f [HobsjZnf 

jj HghsiZn) X {pi}) 

Bh=1^ - ^2 -’ 

n=\ 

28 rrr 1 m2 

Ch = ^ 

n=\ 


{Pi})] 


(54) 

(55) 

(56) 


In above equations subscript obs is refer to observational 
quantities and subscript th is for theoretical one. 


4 Cosmological constraints and data fitting 

As it was mentioned, we have introduced an ansatz as equa¬ 
tion 36, that consist of three free parameters. Where coi indi¬ 
cates present time value of (Ode- For more convenience we 
can suppose Oo = /o = 1 and therefore Eq.(36) is reduced to 
[54,58] 

(ODE{z) = (Oi{l+zf. (57) 

Also, the mean square of relative error functions nor¬ 
mally cause the free parameters plane split in two parts. Peo¬ 
ple usually are interested to the regions which < 1, 

where N denotes the amount of observational data. Whereas 
we use Union — 2 data set for SNela, N for supernovae is 
NsNe = 557, and also for OHD, CMB and BAO, one has 
Nqhd = 28, Ncmb = 1 and Nbao = 1 - Since in this work 
three free parameters are appeared, the space of constraints 
has three dimensions. Thence for better expression, one can 
map figures on two dimensions (in fact it is supposed that, 
the free parameters are independent) and their values will be 
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analyzed. The common regions for best fitting of all con¬ 
straints, play key role in this study. Based on above dis¬ 
cussions we plot a couple of free parameters in Figures 4 
- 6 . In Figure 4 we investigate the constraints on f2„,o in 
0)1 p plate, and also for two constraints SNela and OHD 
minimum points of ore distinguished. In Figures 5 us¬ 
ing best value of Oi, the constraints in f2,„o j3 are obtained, 
in a similar way for best value of j3, the behavior of con¬ 
straints in 0)1 f2„,o surface will be shown. Let us, return our 
attention to figure 4 again. For Q,„q = 0.2, the CMB, BAO 
and OHD have an overlap region, but they are not in agree¬ 
ment with SNela results. Also for a different quantity, the 
SNela and OHD results could be in agreement with together. 
This different behavior of constraints indicates that if one 
wants to compare theoretical results with observations, it is 
better the greatest set of constraints, to be considered. For 
more investigation about overlaps and the effects of indi¬ 
vidual observations, we plot figures 7 and 8 . In figure 7 the 
behaviour of Xj = ZsNe + Xohd + Xcmb + xIao and Xt = 
+ XcMB + Xbao AXj = 3.53,6.25,8.02 are com¬ 
pared. Also in figure 8 to investigate degeneracy one can 

consider Xj = ZsNe + Xohd + Xcmb + Xbao and Xt = Xohd + 
Xcmb +Xbao for ^Xt = 0.1,0.2,0.3. These two figures in¬ 
dicate that although the importance of individual OHD data 
surveying (in comparison SNela, CMB and BAO) is not so 
important, but it decreases the degeneracy between free pa¬ 
rameters of the model. From figures 7 and 8 , it is obvious 
that a collective of four constraints has completely differ¬ 
ent results in comparison to even three constraints. In the 
following, by means of observations, we use some custom 
quantities which are considered for better estimation of the¬ 
oretical parameters of the model. Since all free parameters 
of the model are independent, the total likelihood function 
could be introduced as 


= -XsNe X .ifoHD X .ifcMB X (58) 

therefore the total x^ function could be achieved as 

Xt = XsNe + Xohd + Xcmb+ Xbao- (59) 


It is considerable to attain the maximum amount of the prob¬ 
ability and the minimum value of we should minimize 
Xj- Also it should be noted, in 59 all components have same 
weight. So the likelihood method is equivalent to this fact 
that, for instance all measurements which lead to CMB is 
equal to a supernova explosion!. Afterwards we return to 
this problem. Another quantity which could be used for data 
fitting process is 


X 


Xt 

Ndof 


(60) 


where subscript dof is abbreviation of degree of freedom, 
and Ndof could be defined as the difference between all ob¬ 
servational sources and the amount of free parameters. Let’s 


explain it in more detail, whereas the amounts of all obser¬ 
vations are 557 4-28-1-1-1-1 = 587, and the number of free 
parameters are 4, by considering Hq, therefore Ndof, is equal 
to 583. Also one knows, the acceptable quantity for x^ is 
1.05. For more convenient, we now define the average rela¬ 
tive error functions as follows 


^2 _ ZsNe 

XSNe - AsNe ’ 

.y2 _ XqHT) 

XOHD Nqhd ’ 
2 

y 2 _ )^CMB 

ACME WcMB ’ 

y2 _ Zbao 
/BAO Wbao ■ 


( 61 ) 

(62) 

(63) 

(64) 


Finally we can introduce Xm function, which is equal to 
maximum of X^ functions and it could be considered as 

Xm = max of (XsNe’ Xohd ’ Xcmb >Xbao) ■ (^5) 

In fact the Xm function could be considered as a criterion of 
accuracy for the models. Now we want to compare the be¬ 
havior of Xm ^nd x'^ functions. Without loss the generality 
of the model, one can plot the three dimensional shape of 
Xm und x'^, versus free parameters of the model. These di¬ 
agrams help us to find out the best estimation of the free 
parameters in comparison to observations; for more clar¬ 
ity one can see the figure 9. In this figure, the first dia¬ 
gram shows the minimum of Xm ^^d X^ versus f2„,o. Also 
in two latest diagrams of figure 9, the minimum points are 
drown based on a)i and j3 respectively. By comparison the 
behavior of these relative error functions in Figure 9 one 
can realize that, there are more points (or neighborhood) 
in which X^ < but Xm exceeds 1.05. In fact this behav¬ 
ior was predictable, because in definition of X^, we use the 
contribution of all observational data set. So, for example 
the Xcmb deviation of best fitting results, could be recom¬ 
pense by SNela data abundance. We will return to this draw¬ 
back, after some discussion about likelihood and relative er¬ 
ror functions. For more illustration, we portrait the differ¬ 
ent surfaces of three dimensional, (f 2 jno, , j3), to (Qm 0 ,P), 
(a)i,f2nio) and (a)i,j3) surfaces, which are brought in Fig¬ 
ures 10, 11 and 12. Diagrams B and C are related to (XT)min, 
where the subscript min, shows the minimum value of Xx ■ 
It is notable, in a three dimensional space of free parame¬ 
ters, the confidence levels 68.3%, 90% and 95.4% are pro¬ 
portional to AXx = 3.53, AXx = 6.25 and AXj = 8.02 sur¬ 
faces respectively where AXj = xh~ (XT)min- In diagram 
B of Figures 10, 11 and 12 the counter lines of confidence 
levels are drown and in diagram C, both the Xm surfaces 
and counter lines are brought for more comparison. From 
diagram C it is realized that the confidence level counters 
exceed the Xm regions. From this behaviour it is concluded 
that, the theoretical prediction of CMB shift parameter is 
very greater than it’s observational quantity. As mentioned 
heretofore, when the total mean square error function is in¬ 
troduced the weight of all constraints was identical and this 
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causes some problems. As a matter of fact, the results of 
likelihood’s parameter, equation 59, the effect of CMB shift 
parameter in comparison to the abundant SNela data set is 
ignored. For more information, one can see Table 2 and defi¬ 
nition of Ndof- To overcome these problems, we redefine Xj 
as bellow 

Xt = ZsNe + ZOHD + 3 XcMB + 3 Xbao ■ (66) 

It should be noted, in data fitting and maximization of proba¬ 
bility quantities these two definitions of Xj, i e- equations 59 
and 66 , have not very different. For justifying this claim one 
can compare Tables 2 and 3 which are related to 59 and 66 
respectively. But in figures which related to confidence lev¬ 
els one can observe that the exceeding of confidence lev¬ 
els are reduced, therefore the re-weight of some constraints 
can improve the behaviour of the model. For more clarifica¬ 
tion one can refer to figures 13, 14 and 15. Now by means 
of 66 , we margin the likelihood .if (X 2 nio, , j3) w.r.t (Oi, j3 

and f2mo respectively. Also the relative probability functions 
.if (f2nio, j3), Jf(tUi,f2n,o) and.if(a)i,j3) in two dimensional 
confidence levels 68.3%, 90% and 95.4% are plotted in Fig¬ 
ure 16. For more investigations, we will draw the one di¬ 
mensional marginalized likelihood functions .if(f 2 mo) ver¬ 
sus f2n,o, SX’icoi) based on coi and -Sf(j3) versus j3 in figure 
17. Meanwhile in Table 4 one observes the quantities which 
maximize the marginalized likelihoods using different confi¬ 
dence levels by means of confidence levels ffi = 68.3% and 
(72 =95.4% . 


4.1 Typical example 

Now, we define an effective dark energy as combination of 
dark energy Pde and scalar field density as poE = Pde + 
p<p/f{<p)- So, the Friedmann equation is rewritten as 

3H^ = f{(p) (pm -f Pde^ ■ (67) 

An useful parameter in this study is energy density parame¬ 
ter £2. Here Qoe and respectively will be taken equal to 
^DE = f{(p)PDE/Pc and Qm = f{(p)Pm/Pc, in which Pc is 
the critical energy density which is defined as Pc = 3H^. As 
a result, from the Friedmann equation we have Qqe + = 

1 . 

To obtain energy conservation equations for effective dark 
energy one can achieve the following results 

^(/(^)P£'£) + 3Hf{<p){\ + (ODE)PDE = 7Pmf{<p), (68) 

^(/(^)Pm) +3///((p)(l+7)Pm = -7Pm/(9), (69) 

SO that the effective pressure of dark energy is defined as 
Pde = Pa + P<p/f{(p), and one has the effective dark energy 
equation of state parameter as code = Pde/Pde- Also 7 is 


the matter equation of state parameter which is defined as 
7 = Pm/Pm- For 7 = constant, integrating of Eq.(69) results 
in the following relation for cold dark matter energy density 


Pm — 


P** 

rem 


(70) 


a3{i+r)y(i+y)(^)’ 

where p®,, = fQ^^\<p)Pm- In this step, we suppose that the 
effective dark energy could be defined as ADE, in other 
word we assume that 


3n^ 

Pde = Pade = yj, 


(71) 


where n is a numerical constant and T is cosmic time and 
therefore Qde is obtained as Qde = f{(p)n^/H^T^- Taking 
this assumption and using Eq.( 68 ), the equation of state pa¬ 
rameter of effective dark energy could be acquired as 


CODE = — 1 + 


2 1 Qde , f{(p) 




(72) 


3n^ficp) 3Hf{cp) 

where r is ratio of cold dark matter and effective dark en¬ 
ergy, namely r = Pm/PDE = ^m/^DE- The interaction term 
in this model generates an extra term for code, which can 
justify the phantom divide line crossing. By definition an 
ansatz for cOgA, it can be considered as 

cOgA 3-i = coq C0i{lz)^■ (73) 

Eor fitting the free parameters for ADE in an external scalar 
field interaction model, we use the 557 Union-2 sample database 
of SNela, and p^j = Pradiation “t“ Pbaryon % Pdarkmatrer- There¬ 
fore in this case the Eriemann equation is as 


3H^ 


-■fi(p)[pm + pDE)- 


Combining Eqs.(68)-(71), give 


3H^ = 


V( 


p^ 

rem 


3n^\ 


(74) 


(75) 


a3(i+y)y{i+r)(^) pZ 7 ’ 

where Pg,„ is the effective energy density of matter at the 
present time. Whereas the 557 Union-2 sample database have 
collected from red shift parameter to various SNela, there¬ 
fore we rewrite E = H/Hq versus z as 


E^ = 


fo(l +z)^ + (1 +z)3““exp{3^[(l - 1 ]} 


.(76) 

ro +1 

To achieve the best fit for free parameters based on subsec¬ 
tion 3.1a minimization method leads to 


ZL„(^ = 1-1; ®i =-1.65; j3 =-2.25), 
2 

Xiin = ^ - TT = 542.75, 


C 


B 


Po = -^ =43.1089. 


(77) 

(78) 

(79) 


where implies Xm/dof= x]n„Jdof = Q.9%\{dof = 553). 
In figure 1, we show a comparison between theoretical dis¬ 
tance modulus and observed distance modulus of supernovae 












Fig. 1 The observed distance modulus of supemovae (points) and the 
theoretical predicted distance modulus (red-solid line) in the context of 
ADE model. 



Fig. 3 The plot shows the evolution of effective dark energy parame¬ 
ter, CODE, versus z, for co^) = 1.1, £ 0 i = —1.68 and j8 = —2.25. 



Fig. 2 Contour plots for the free parameters (Oi and /3 , shows that the 
best value for these parameters are — 1.86 < tai < — 1.62 and —2.27 < 
P < -0.73. 


here /o is the constant of integration. Whereas 


f{(p) 

f{<P) 




A significant result of observational data is accelerated ex¬ 
pansion of the Universe. A good cosmological model should 
be able to describe this acceleration. An useful quantity to 
investigate this property of the Universe is deceleration pa¬ 
rameter which defined as g' = —1 Using Eqs.(67), 

(68) and (69), one achieves the deceleration parameter gives 


^ ^ 1 ^ — Oi(l-t-z)^ 

( _ Do _ \ 


(81) 


where Dq is the constant of integration. It is clearly seen 
that for 0)0 = 1.1, 0 ) 1 =—1.68, j3 =—2.25, (which have 

obtained from data fitting processes) q <0. 


data. The red-solid line indicates the theoretical value of dis¬ 
tance modulus, jith, for the best value of free parameters 
0)1 = —1.65, 0)0 = 1.1 and j3 = —2.25. This shows that the 
model is clearly consistent with the data since /dof = 1. 
Fig.2 show contour plots for the free parameters 0)i and 
j3, it is shown that the best value for these parameters are 
— 1.86 < 0)i < —1.62 and —2.27 < < —0.73 in which for 

stability condition > 0, we have taken the interface be¬ 
tween green and yellow sector, 0)i = —1.68. 


The evolution of effective dark energy parameter, (Ode, ver¬ 
sus z, for 0)0 = 1.1, 0)1 = —1.68 and j3 = —2.25 have been 
shown in Fig.3. This show that by growing z the parameter 
get into the phantom phase. 

Here 0)o, 0)i and j3 are free parameters of the model which 
obtained from data fitting. It is clear that, if the form of dark 
energy density is given the coupling function, f{<p), could 
be easily determined. For instance by using Fqs.(68), (69) 
and (7 1) one can obtain 


/((p)=/ofV3“Oexp 



13+2 J’ 


(80) 


5 Conclusion and Discussion 

Interacting models which contain an external interaction be¬ 
tween matter and scalar fields attract more attentions. Such 
mechanisms are capable to suppress the fifth force and also 
are in good agrement with observations. Using such pow¬ 
erful mechanism we have attained some cosmological pa¬ 
rameters consist of coincidence and deceleration parame¬ 
ters. For instance based on Table 2, and equations 32 and 
34 it is clear that r(z) is a decreasing function and q has 
taken negative values for different amounts of z. Consid¬ 
ering a suitable ansatz for FoS parameter of effective dark 
energy dimensionless Hubble parameter is obtained. So by 
means of SNela, CMB, BAO and OHD data sets the min¬ 
imum value of X^ for the free parameters of the model are 
achieved. To estimate the free parameters of an ansatz for the 
effective dark energy equation of state, the whole of the ob¬ 
servational data sets have been considered. For more details 
one can compare the results of figures 7 and 8, and the re¬ 
sults of typical example, subsection 4.1. Also for getting bet¬ 
ter overlap between the counters with the constraint Xm — 3, 
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the Xj function have been re-weighted. Meanwhile the rela¬ 
tive probability functions have plotted for marginalized like¬ 
lihood (f2nio, Wi, j3) according to two dimensional confi¬ 
dence levels 68.3%, 90% and 95.4%. In addition the value 
of free parameters which maximize the marginalized likeli¬ 
hoods using above confidence levels have obtained. Based 
on above discussions a couple of free parameters in figures 
4-6, have plotted. In figure 4, the constraints on f2„,o in 
Oi p plate have investigated; and also, for two constraints 
SNela and OHD minimum points of have distinguished. 
In figures 5, using best value of Oi, the constraints in Q^q j3 
surface are obtained. In a similar way, for best value of j3, 
the behavior of constraints in (0i Q„,o plane have shown. 
Also based on figure 4, for Q,no = 0.2, the CMB, BAO and 
OHD have an overlap region, but they are not in agreement 
with SNela results; where one possible explanation would 
be an incompatibility among the data sets. Also, for a dif¬ 
ferent values one can find a region which SNela and OHD 
are in more agreement against CMB and BAO. This differ¬ 
ent behavior of constraints indicates that, if one wants to 
compare theoretical and observational results, it maybe bet¬ 
ter the greatest set of constraints to be considered. For more 
investigation about overlaps and the effects on individual ob¬ 
servations, the figures 7 and 8 , have been plotted. In figure 
7, the behaviour of Xj = ZsNe + Xohd + Xcmb + xIao 
Xt = XsNe + XcMB +Xbao AXj = 3.53,6.25,8.02, have 
compared. Also in figure 8, we have considered Xj = XsNe + 

XoHD + Xcmb + Xbao Xj = Xohd + Xcmb + Xbao’ 
AXt = 0.1,0.2,0.3, to investigate the degeneracy. These two 
figures indicate that although the importance of individual 
OHD data surveying in cosmological investigations (in com¬ 
parison SNela, CMB and BAO) is not so important but it 
decreases degeneracy between free parameters. 

Acknowledgements H. Sheikhahmadi would like to thank Iran’s Na¬ 
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Table 1 In this table from left to right z, H{z){kms^^ Mpc^^), it’s un¬ 
certainty Gif {km Mpc^^) in measurement and related references 
(by considering the technique which is used) are collected, respec¬ 
tively. 



H{z) 

Oh 

References 

Techneque 



0.070 

69 

19.6 

[65] 

SDSS DR7; 0 < z 

<0.4 

0.100 

69 

12 

[66] 

ATC;0.1 <z< 

1.8 

0.120 

68.6 

26.2 

[65] 

SDSS DR7; 0 < z 

<0.4 

0.170 

83 

8 

[66] 

ATC;0.1 <z< 

1.8 

0.179 

75 

4 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.199 

75 

5 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.200 

72.9 

29.6 

[65] 

SDSS DR7; 0 < z 

<0.4 

0.270 

77 

14 

[66] 

ATC;0.1 <z< 

1.8 

0.280 

88.8 

36.6 

[65] 

SDSS DR7; 0 < z 

<0.4 

0.350 

163 

5.6 

[68] 

SDSS DR7 LRGs; 

z=0.35 

0.352 

83 

14 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.400 

95 

17 

[66] 

ATC;0.1 <z< 

1.8 

0.440 

82.6 

7.8 

[69] 

WiggleZ-tH(z); z 

< 1.0 

0.480 

97 

62 

[70] 

CMB-tOHD; 0.2 < 

z< 1.0 

0.593 

104 

13 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.600 

87.9 

6.1 

[69] 

WiggleZ-tH(z); z 

< 1.0 

0.680 

92 

8 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.730 

97.3 

7.0 

[69] 

WiggleZ-tH(z); z 

< 1.0 

0.781 

105 

12 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.875 

125 

17 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

0.880 

90 

40 

[70] 

CMB-l-OHD; 0.2 < 

z < 1.0 

0.900 

117 

23 

[66] 

ATC;0.1 <z< 

1.8 

1.037 

154 

20 

[67] 

OHD-tCMB; 0 < z 

< 1.75 

1.300 

168 

17 

[66] 

ATC;0.1 <z< 

1.8 

1.430 

ill 

18 

[66] 

ATC;0.1 <z< 

1.8 

1.530 

140 

14 

[66] 

ATC;0.1 <z< 

1.8 

1.750 

202 

40 

[66] 

ATC;0.1 <z< 

1.8 

2.300 

224 

8 

[63] 

BAD; 0.7 < z < 

2.3 


3. 


2 , 


1 
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Fig. 4 The counter lines of XsNe = 557 (brown), Xohd ~ (green), 
XcMB ~ * (red), and Xbao ~ 1 ■0{blue) for of Qi„q = 0.26 are plotted. 
Also for two constraints SNela and OHD minimum points of x^ are 
distinguished. The dashed lines refer to the counter lines which are 
greater of the minimum points only unity. 
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Fig. 5 The counter lines of = 557 (brown), Xohd ~ (green), 
XcMB =1-0 (red), and Xbao ~ \-Q{blue) for tt)i = —1.1 are plotted. 
Also for two constraints SNela and OHD minimum point of x^ are 
distinguished. The dashed lines refer to the counter lines which are 
greater of the minimum points only unity. 



<Ji 

Fig. 6 The counter lines of XsHe ~ 557 (brown), Xohd ~ (gr^eri), 
XcMB =1-0 (red), and Xbao ~ \^-^{blue) for P = —0.25 are plotted. 
Also for two constraints SNela and OHD minimum point of x^ are 
distinguished. The dashed lines refer to the counter lines which are 
greater of the minimum points only unity. 


Table 2 In the table, the quantities related to minimum point of Xj ~ 
2;sNe + 2 (ohd + 2;cmb + ^BAO introduced. 
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Zbao 
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16 X 10-^, 
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(Zt)™!! 
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Fig. 7 In this figure the behaviour of Zt = 2;sNe + Zohd + ^^cmb + 
Zlf^o (dashed counters)and Xr = + ^cmb + ^bao (^o^d counters) 

for AXy ~ 3.53{innerloops),6.25(middleloops),S.02{outerloop) 
are compared. The minimum points of these two Xj functions are dis¬ 
tinguished by Solid points. 



Fig. 8 In this plot we consider Xj = ZsNe + Xohd + Xcmb + xIao 
(dashed counters) and Xj = Xohd + ^cmb + ^bao (^‘’''d counters) for 
AXj = Q. {{innerloops), 0.2{middleloops), 0.3{outerloop) to investi¬ 
gate degeneracy in this work. This Figure and Figure 7 indicate that 
although the importance of individual OHD data surveying in cosmo¬ 
logical investigations (in comparison SNe la, CMB and BAO) is not 
so important but it causes decreasing degeneracy between free param¬ 
eters of the model. The minimum points of these two Xt functions are 
distinguished by Solid points. 
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Fig. 9 In above diagrams the minimum quantity of X^ (blue line) and 
Xm (red-upper-line) versus tt>i and /3 parameters have been drown 
respectively. 


Table 3 This table is related to minimum point of Xj = /(sNe + ^ohd + 
32;cmb + ^Xbaq- 




COi 

■^mO 

i^^BAO 

-0.239 

-1.051 

0.272 
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Fig. 10 In diagram (A), the image of Xm < 1 on the 
(f2„,o,/3) surface are portrait. In (B) the minimum point 
of Xt = ZsNe + XoHD + XcMB + xLo ^nd the shadow of 
AXj ~ 3.53(innerloop),6.25{middleloop),9i.02{outerloop) .surfaces 
on the (f2„,o,/3) plate, are plotted. In part (C) both diagrams (A) and 
(B) are brought to compare the results. 



Fig. 12 In diagram (A), the image of Xm < 1 on the 
surface are portrait. In (B) the minimum point 
of Xf = XsNt + XoHD + 2;cmb + ^BAO ^nd the shadow of 
AXj = 3.53{innerloop),6.25(middleloop),S.02{outerloop) surfaces 
on the (cOi.jS) plate , as counter lines, are drawn. In part (C) both 
diagrams (A) and (B) are compared. 



Fig. 11 In diagram (A), the image of Xm < 1 on the 
(C()i,f2,„o) surface are portrait. In (B) the minimum point 
of Xt = + ^OHD + 2;cMB + xIao ^nd the shadow of 

AXj ~ 3.53{innerloop),6.25(middleloop),?i.02{outerloop) surfaces 
on the {a)i,Q„,o) plate, are drawn. In part (C) both diagrams (A) and 
(B) are considered for more comparison. 


Fig. 13 In diagram (A), the image of Xm < 1 on the {Q,„o,l3) 
surface are portrait. In (B) the minimum point of Xj ~ 
^SNe + ^OHD + 3 ^cmb + '^xIaq and the shadow of AXt = 
3.53(innerloop), 6.25(middle loop), 8.02(outerloop) surfaces on the 
(i2„,o, P) surface, as counter lines, are drawn. In part (C) both diagrams 
(A) and (B) have been brought for more comparison. 
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Fig. 14 In diagram (A), the image of Xm < 1 on the (C0i,f2„,o) 
surface are portrait. In (B) the minimum point of Xt ~ 
:i:sNe + XoHD + 3;i;cmb + and the shadow of Axk = 

3.53{innerloop),6.25(middleloop),H.02{outerloop) surfaces on the 
(cOi, f2„,o) plate, as counter lines, are drawn. In part (C) both diagrams 
(A) and (B) are brought for comparison. 



Fig. 16 In the above two dimensional likelihood diagrams, the 68.3% 
confidence level (dotted line), 90.0% confidence level (green-dashed- 
line) and 95.45% confidence level (red-solid-line) after marginalization 
on the coi, j8 and Q,„q free parameters are plotted. Note that in this fig¬ 
ure we use Xt = ZsNe + Xom + ^Xcmb + "ixlAO and also the shadow 
of AXj ~ 3.53{innerloop),6.25{middleloop),S.02{outerloop) sur¬ 
faces. 



Table 4 In this table the quantities which maximize the relative prob¬ 
ability functions .if(£0i) and J5f(j8) using confidence levels 

Cl = 68.3% and CT 2 = 95.4% are calculated. The data sets are includes 
of SNela, CMB, BAO and OHD in which the weight of Xcmb and Xbao 
in Xrotal function is the coefficient 3. 
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0.44 
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Fig. 15 In diagram (A), the image of < 1 on the (coi.jS) 
surface are portrait. In (B) the minimum point of Xt ~ 
:S:SNe + ^OHD + ^XcMB + ^xIao and the shadow of AXj ^ 
3.53{innerloop),6.25{middIeloop),S.02[outerIoop) surfaces on the 
(cOijjS) plate, as counter lines, are drawn. In part (C) both diagrams 
(A) and (B) are collected for comparison. 
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Fig. 17 In above diagrams the relative likelihoods are plotted. The 
68.3% and 95.4% confidence levels are distinguished as brown dashed 
line and red dot-dash line respectively. It should be noted that in 
these figure we use Zj ^ Zsnc + ^ohd + the 

best fits of the free parameters are as j3 = —0-24 ^q23’^oS’ ~ 

^.+0.08,+0.156 _ p. ^^^+0.01,+0.02 

^■^^-0.08,-0.16 “ ^■"^^"^-0.01,-0.02- 















































